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Introduction 


M.  I.  Freidlin  ([9],  [10])  has  recently  introduced  probabilistic 
techniques  Into  the  study  "wavefront  propagation"  for  systems  of  reaction- 
diffusion  PDE.  The  motivating  Idea  is  that  should  a  reaction-diffusion  system 
possess  only  a  single  unstable  and  a  single  stable  equilibrium,  then  the 
solution  u  of  the  system  will  presumably  tend  to  "switch"  for  large  times 
from  near  the  former  to  near  the  latter  state.  A  mathematical  problem  Is  then 
to  describe  this  transition,  ideally  in  terms  of  simpler  quantities  than  those 
governing  the  full,  detailed  behavior  of  the  entire  system  of  PDE.  More 
precisely,  to  study  the  reaction-diffusion  problem  for  large  times  of  order 
e  \  Freidlin  suggests  a  e  *  rescaling  in  the  space  and  time  variables;  so 

e 

that  our  attention  turns  to  the  solutions  u  of  certain  c-dependent  systems 

g 

of  PDE.  We  then  hope  to  show  that  as  e  — »  0,  the  functions  u  converge  in 
some  region  G  c  Rn  x  [0,oo)  to  the  stable  point,  and  in  the  opposite  region 
[IRn  x  [0,»)]\G  to  the  unstable  point.  We  simultaneously  hope  to  describe 
geometrically  or  analytically  this  set  G,  whose  boundary  we  envision  as  a 
spreading  wavefront  separating  regions  with  quite  different  limiting  behavior. 

This  paper,  which  is  an  extension  to  systems  of  earlier  work  [6]  on 
single  equations,  brings  to  bear  purely  PDE  techniques  to  this  problem, 
especially  the  theory  of  viscosity  solutions  on  Hamilton- Jacobi  equations, 
due  to  Crandall -Lions  [3].  The  connection  with  the  foregoing  discussion  is 
that,  the  region  G  alluded  to  above  is  the  set  where  the  solution  J  of  a 
certain  Hamilton- Jacobi  equation  is  negative.  Our  procedure  for  understanding 
the  limiting  behavior  of  the  solution  of  the  reaction-diffusion  system  of  PDE 
is  thus  first  of  all  to  build  an  appropriate  Hamiltonian  H  out  of  the  data 
given  In  the  problem,  second  to  solve  the  resulting  Hamilton- Jacobi  equation 
for  J,  and  last  to  demonstrate  the  different  limiting  behavior  of  the 
solutions  uc  of  the  scaled  system  on  the  sets  {J<0>  and  {J>0}.  We 
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Informally  regard  the  Hamiltonian  as  controlling  somehow  the  rate  of  instabi¬ 
lity  of  the  unstable  point.  We  are  thus  able  to  characterize  the  asymptotic 
behavior  of  the  "complicated"  reaction-diffusion  system  in  terms  of  the 
"simple"  Hamilton-Jacobi  equation.  This  possibility,  first  Identified  by 
Freidlin  [9]  in  rather  different  terras,  is  attractive,  but  of  course  requires 
for  this  implementation  many  structural  assumptions  on  the  nonlinearities, 
which  we  list  below.  It  would  of  course  be  quite  interesting  to  extend  our 
results,  or  at  least  the  point  of  view  espoused  above,  to  systems  with  more 
general  nonlinearities. 

More  precisely  now,  we  intend  to  investigate  the  scaled  reaction- 
diffusion  system: 

u£  t  =  edkAu£  +  “  fk(uG)  in  Rn  x  (0,co) 

(1.1), 

£ 

u£  =  gk  on  Rn  x  {0}  (k=l,  • • • ,m) . 

Here  the  constants  d^  (l£k£m),  and  the  functions 

_.n  _m  ,  „  _m  _m 
g:  R  — »  R  and  f:  R  — »  R  , 

are  given,  where  we  write  g  =  (g, ),  f  =  (f,,***,f  ).  The  unknown  is 

l  m  l  m 

u  ■  (u, , •• *,u  ).  We  will  assume  that 
i  m 

d^  >  0  (k=l, •••,m), 

4 

and  that  the  functions  g,  f  ere  smooth,  bounded  and  Llpschitz. 

In  addition  we  suppose  that 

(1.2)  gk  *  0  (k=l, • • • , m) 

and 

(1.3)  GQ  ■  (gk  >  0)  (k=l , • • • , m) 
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Is  a  bounded,  smooth  subset  of  Fn.  Under  these  assumptions  there  exists  a 

£ 

unique  smooth  solution  u  of  the  PDE  (1.1)  ,  with 

c 

u£  >  0  in  Fn  x  (0,»)  (k=l, • • • , m) . 

Our  essential  assumptions  all  concern  the  reaction  term  f.  First  of 
all  we  suppose 

(FI)  f(0)  =  0  ; 

and  alsc 

'  ,kl"i''"’Vr0'Vi'"',“il  >  0  lf 

(F2) 

Uj,  ♦••,Ujc_j,Ujc+j,  •••,u  t  0  and  u^  >  0  for  some  index  l*k. 

Consequently  the  vector  vleld  f  points  strictly  inward  along  the  boundary 
of  the  positivity  set 

II  *  (u  €  Rm|u, >0,  •  •  • ,  u  >0} , 

1  i  m 

except  at  the  point  0,  which  is  an  equilibrium  point  for  the  system  (1.1)  . 

£ 

To  ensure  that  our  solutions  u  do  not  become  unbounded  as  c  — »  0,  we 
further  hypothesize  that 

there  exists  a  constant  A  such  that 

fk(u)  S  0  (k=l, • •*,m). 

if  u  e  II  and  u^  £  A. 


C  ■  Df ( 0 ) , 


Df  denoting  the  gradient  of  f.  Thus 


ckl  *  'k.u/01 


We  assume 


(ISk, ISm). 


Ckl  >  0 


(ISk, ISm), 


fk(u)  *  cklul 


(u«n,  k«l,»»»,m), 


where  we  employ  the  standard  summation  convention. 
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for  systems  of  the  Kolmogorov-Petrovski i-Piskunov  nonlinearity,  discussed  In 

[6] . 

Our  main  result,  Theorem  1,  asserts  that  under  hypotheses  (FI)  -  (F5) 

g 

u  (x,t)  converges  as  e  — *  0  to  zero  or  not  depending  on  whether  J(x,t)  >  0 
or  J(x,t)  <  0,  the  function  J  satisfying  a  Hamilton- Jacobi  PDE  whose 
structure  we  now  describe.  Given  p  €  Fn,  define  the  m  x  m  matrix 

B(p)  *  dlag(***,  dk|p|2,***) 

and  then  set 

(1.4)  A(p)  ■  B(p)  +  C. 


Now  the  matrix  A(p)  has  positive  entries,  and  so  Perron-Frobenius  theory 
asserts  that  A(p)  possesses  a  simple,  real  eigenvalue  A°  =  A°(A(p)) 
satisfying 

Re  A  <  A0 


for  all  other  eigenvalues  A  of  A(p).  Let  us  define  then  the  Hamiltonian 

(1.5)  H(p)  «  A°(A(p))  (P€Rn). 

See,  for  instance,  [7]  for  a  review  of  the  various  properties  of  H,  and  in 
particular  a  proof  that  H  is  convex.  We  additionally  set 

(1.6)  L(q)  *  sup  (p»q-H(p)); 

pcR 

L  is  the  lAgranglan  associated  with  H.  Finally  we  define  for  each  point 
(x.t)  e  Rn  x  (0, * )  the  action  function 


(1.7) 


J(x,t)  *  inf  J  L(z(s)  )ds|z(0)  €  GQ,z(t)  =  xj-, 


the  infimum  taken  over  all  absolutely  continuous  functions 


satisfying  the  stated  initial  and  terminal  conditions.  As  we  will  see,  J 
turns  out  to  be  the  (unique)  solution  of  the  Hamllton~Jacobl  equation 


+  H(DJ)  =  0 

in 

Rn  x  (0,®) 

(1.8) 

J 

=  0 

on 

o° 

X 

o 

J 

k 

-  +  ® 

on 

lnt(Rn-GQ) 

in  the  viscosity  sense  (cf.  Crandal 1 -Lions  [3],  Crandall -Evans-Lions  [2], 
etc. ) . 


Theorem  1.  Under  hypotheses  (Fl)-(FS)  we  have 
(1.9)  '*-  c 

and 


11m  u  -  0 
e-X) 


uniformly  on  compact  subset  of  {J>0} 


(1.10)  lim  inf  u  >  0  (k«l,*»*,m)  uniformly  on  compact  subset  of  {J<0> 

e-X) 


We  loosely  interpret  this  theorem  as  describing  how  the  Hamiltonian  H, 

which  depends  upon  both  C  *  Df(0)  and  the  diffusion  constants  d  .••♦,d 

1  m 

controls  the  instability  of  the  equilibrium  point  u  =  0. 


Remark  We  should  note  also  that  it  is  possible  to  refine  conclusion  (1.10) 
by  making  further  assumptions  on  the  behavior  of  the  vector  field  f  in  IT. 
As  in  Freldlin  [9],  we  may  for  instance  assume  that  there  exists  a  unique 
equilibrium  point  a  €  II  which  is  asymptotically  stable  for  the  flow 
generated  by  the  vector  field  f,  as  in  the  following  illustration. 
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Under  various  fairly  stringent  technical  assumptions,  it  is  then  possible  to 
show  that 

(1.11)  11m  uE  =  a  uniformly  on  compact  subsets  of  {J<0>. 

e-»0 

We  may  consequently  think  of  the  boundary  of  the  set  G  a  {J<0}  as  a 
"propagating  wavefront"  in  the  sense  explained  above. 

Assertions  (1.9),  (1.11)  are  attractive,  but  the  sufficient  conditions 
we  know  to  improve  (1.10)  to  (1.11)  are  rather  unsatisfactory  technically. 

(For  instance,  we  could  assume  that  the  diffusions  constants  are  all  equal  and 
that  f  has  a  convex  Lyapunov  function  associated  with  the  stable  point  a, 
etc.).  For  this  reason,  we  will  not  here  formulate  any  precise  assertions 
leading  to  (1.11),  but  instead  refer  the  reader  to  Freidlin  [9],  [10]  for 
more  information  on  this  point.  □ 

Our  paper  Is  organized  so  that  the  proof  of  Theorem  1  appears  in  £3,  £2 
being  devoted  to  some  preliminary  estimates.  The  main  idea,  following  [8], 
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is  to  perform  a  change  of  variables,  after  which  we  send  c  — »  0.  This  approach 
is  greatly  complicated  for  the  present  case  of  systems,  since  we  lose  the 
maximum  principle  and  consequently  many  of  the  estimates  available  in  [6]  for 
the  scalar  case.  We  depend  Instead  upon  some  recent  techniques  introduced  by 
Ishii  [11]  and  Barles-Perthame  [1]  in  Hamilton- Jacobi  theory,  which  allow  us 
to  proceed  with  only  supremum-norm  bounds.  This  is  a  great  advantage  since 
better  estimates  seem  to  be  unobtainable,  but  the  price  is  a  far  greater 
intricacy  in  some  of  the  convergence  proofs.  We  will  encounter  for  instance 
troubles  in  interpreting  in  Just  what  sense  certain  solutions  take  on  their 
initial  values.  An  appendix  (f4)  develops  some  useful  theory  on  such 
questions,  which  we  will  need  in  f3. 

Finally  we  let  us  note  that  although  Freidlin’s  work  has  greatly  inspired 
us,  we  believe  that  the  PDE  techniques  developed  here  (and  in  [6])  provide 
information  which  is  not  at  all  clear  from  the  probabilistic  viewpoint.  We 
for  instance  do  not  require  the  fairly  specific  structural  assumptions  on  the 
nonlinearity  f  utilized  in  [9,  p.467]. 
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2.  Preliminary  estimates 


Henceforth  we  always  suppose  hypotheses  (F1MF5)  to  be  In  effect 


Lemma  2. 1  There  exists  a  constant  such  that 

(2.1)  0  <  u£  £  Cj  __-ln  Fn  x  (0.®) 


for  k=l. 


,  m  and  each  c 


>? 


Proof .  Choose  a  smooth,  bounded,  Lipschitz  vector  field 

t:  Rm  -*  R™ 


such  that 

(2.2) 


?(u)  =  f(u) 

•  and 

?k(u)  i  0 


(u  e  n) 

IP  <  0,  k®l,***,m; 


this  is  possible  in  light  of  (F2).  Now  let  u  be  the  unique,  smooth 
solution  of  the  system 

V“C)  In  Rn  x  (0,») 

*  gk  on  Fn  x  {0>. 

Choose 


7):  F  — »  R 

to  be  smooth  and  convex,  with 


(2.4) 


’  V  =  0 
.  tj  >  0 


on  (0,«*) 
on  (-o#,0) 


Then  for  each  k*l .  a, 
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at  /  n’’('£]dt  =  ' c  J  Ki2dx 

R  R 


Since 


S  0. 


by  (F3) . 


h 


Trtgk)dx  =  o, 


we  discover 


\ 


G  £  8  in  Rn  x  [0,«). 


3.  Observe  now  that  hypotheses  (FI),  (F2)  and  (F4)  imply  f^tu)  ^  0  if 
£  0  is  small  enough  (k=l,**»,m).  Thus  the  strong  maximum  principle  for 
parabolic  equations  implies 


if  >  0  in  IRn  x  to,®),  (k=l,***,m)  o 


We  next  introduce,  following  [5],  16],  [7],  [8],  etc.  the  new  funct ions 


(2.5) 


=  -  clog  u£  in  Rn  x  (0,«)  (k=l,***,m). 


C  G  G 

A  calculation  shows  that  v  =  (v,,***,v  )  satisfies  the  system 

1  m 


(2.6) 


vk,t  -  cdkivk  *  dklDvklz 


rk(u'>  n 

in  Rx(0,») 


c 


-  c  log  gk  on  Gg  X  {0> 


V.  =  +  09 

k 


on  int(Rn-GQ)  x  {0>. 


Note  also  that  in  view  of  Lemma  2. 1  we  have 


(2.7) 


v  *  -  c  log  C. 


In  Rn  x  [0,®) 
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for  k=l , 


,m  and  each  e  >  0. 


Lemma  2 . 2  For  each  compact  subset  Q  c  c  (Rnx(0,»))  v  (Gq  x  [ 0 , ») )  there 
exists  a  constant  ( Q)  such  that 

|vj|  £  C2(Q)  in  Q 

for  k  =  l,***,m  and  each  e  >  0. 


Proof  1.  Let  0  denote  any  open  subset  of  Rn  x  (0,»).  Suppose  that  0  is 
a  smooth  function  satisfying 


(2.8) 


*t  ~  edkA*  ♦  dk  |D0|2  *  A 


in  0, 


where  A  >  0  is  a  positive  constant  to  be  selected  later.  Assume  now  that 
£ 

vk~0  has  a  maximum  at  some  point  (x^.t^)  €  0.  Then  at  the  point  (x^.t^)  we 
have 


But 


° s  (v*]t  -  «v(vH 

f  (uc) 

<  -  — - A  by  (2.6)  and  (2.8). 

c  e 

\ 


Vu  )  s  -  V--V1-0- 


[fk(’**,uk-r0,  Vr”)'V,,,’V1’VVr*')] 


£ 

S  C  uk  ,  according  to  Lemma  2. 1  and  (F2). 
Substituting  above  we  discover 

0  i  C  -  1, 
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a  contradiction  for  X  large  enough. 


2.  We  may  assume  that  for  some  r  >  0  the  ball  B(0,r)  lies  in  CQ.  Set 
then 


♦'  *  T77 2 4  st  4 » 

r  -|x| 


(|x|<r,  t*0) 


P  and  y  picked  as  indicated  below  (cf.  [51,  [7]).  We  compute 


1,2 


4>[  -  edjA*1  +  dk}D^1| 


=  B  -  edL 


2n 


_§M 


[r2-|x|2  (r2-jx|2)3  J 


+  d. 


k/  2  |  i2\4 
(r  I x |  ) 


£  X  in  B(0,  r)  x  (0,») 


(k=l , 


provided  p  >  0  is  large  enough.  If  we  then  adjust  y  to  be  so  large  that 


l  c 
*  2  vk 


on  B(0,r)  x  {0}, 


we  deduce  from  Step  1  that 


e  .  ,1 

vk  S  * 


In  B(0,r)  x  (0,«0. 


Consequently 

(2.9) 


\sc 


in  B(0,r/2)  x  (0,T). 


for  k=l,***,m  and  each  e  >  0,  T  >  0. 


3.  Now  write 


*2  .  iliif  .  gt  .  r 


( |x|>|,t>0). 


We  calculate 


*2  -  cdkA^2  +  dk|D^2J2 


«l*l  T 

t 


2  2n  +  d 


p  -  c  d. 


k  t2 
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iA  in  (Rn-B(0, r/2) )  x  (0,«), 

provided  a,  0  >  0  and  large  enough.  We  now  pick  y  large  enough  to  ensure 

<PZ  *  v£  on  PB(0, r/2)  x  (0,T). 

This  is  possible  owing  to  estimate  (2.9).  Consequently  Step  1  implies 

\  S  *2  in  (Rn-B(0, r/2) )x(0,T) 

for  k=l,*»*,m  and  each  e  >  0,  T  >  0. 

This  bound  and  (2.9)  lead  at  once  to  the  estimate  stated  in  the 
Lemma.  o 
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3.  Proof  of  Theorem  1. 

For  each  (x,t)  e  Rn  x  (0,«)  we  set 

(3.1)  vk,(x,t)  *  11m  Inf  v£(y,s)  (k=l,***,m) 

c-)0 

(y,s)-»(x,  t) 

and 

*  G 

(3.2)  v^fx.t)  h  lim  sup  vk(y,s)  (k*l,**»,m). 

e-)0 

(y,s)+(x,t) 

We  additionally  write 

*  * 

(3.3)  v„  a  min  v  v  *  max  v.  . 

lSk<m  ISkSm 

( 

We  intend  to  show  that  v,  is  a  lower  semi continuous  supersolution,  and  v 
is  a  upper  semi continuous  subsolution,  of  an  appropriate  variational 
inequality  involving  the  Hamiltonian  H  defined  in  fl. 

Proposition  3. 1.  We  have 
(i) 

(3.4)  min(vn  +  H(Dv.),v.)  20  in  Rn  x  (0,«o) 

in  the  viscosity  sense,  and 
(li) 

0  on  Gq  x  (0} 

(3.5)  v,  =  - 

+«  on  (Rn-Gn)  x  {0>. 

Proof  1.  Because  of  (2.7)  it  is  clear  that 

v.  2  0. 

To  demonstrate  that 

v.t  +  H(Dv.)  20  in  Rn  x  (0,eo) 
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in  the  viscosity  sense,  we  fix  a  smooth  test  function  0  and  assume 


(3.6) 


v,-#  has  a  strict  local  minimum  at  some  point 
.  ( xQ >  1 0 )  €  Rn  X  (0,«). 


We  must  prove 

(3.7)  <pt  +  H(D0)  *  0  at  (x0,tQ). 


2.  Let  (0j,***,0k)  be  a  positive  eigenvector  of  the  matrix  Ad^x^tg)) 
corresponding  to  the  principal  eigenvalue  A^(A(D0(Xq,  t^) )  *  H(D0(xQ, t^) ) .  . 
Note  then 


(3.8) 


v.(x,t)  =  lim  inf  min 

c-O  1515m 

(y,  s)->(x,  t) 


(y,s)  +  c  log 


4 


Combining  (3.6)  and  (3.8),  we  deduce  that  there  exists  an  index  k  e  {l,***,m>, 

C  G 

a  sequence  — »  0,  and  points  (x  r,t  r)  such  that 


(3.9)  v, 


;r(*Cr.tCr] 


+  e  log  0  =  min 

r  151 5m 


[v'r[xCMc']  ♦  cr  log  ->  V. 


'VV’ 


(3. 10) 


(vkr  *  cr  l0*  '"J  - 


0  has  a  local  minimum  at 


(xcMCrj. 


Applying  then  the  maximum  principle,  we  obtain  from  (2.6)c  the  estimate 


2  fJuei 

0  S  *t  -  erdkA*  ♦  dk|D^r  +  I 

V 


(3. 12) 


*  4>t  ~  erdkA*  +  dkN|  +  ckJ  exp 


fvk,'"vi1 

+  c.  .  exo  -1 - *- 


at  the  point  jxe,\terj,  the  second  inequality  being  a  consequence  of 
hypothesis  (FS).  Using  (3.9)  and  (3.11),  we  simplify  (3.12)  to  read 

2  *i 

°  S  4>t  +  dk|D0|  +  ckl  jp  +  o(1)  33  er  0 

S  ¥»k 


at  (x0,tQ).  Since 


dklD*l  \  *  ckl  *1  *  (*<Wx0.t0))«)k  -  HW<x0,t0Wk, 


we  deduce  (3.7)  upon  letting  »  0. 


3.  We  next  verify  assertion  (ii).  To  do  so,  we  fix  p  >  0  and  select 
<  €  c“(IRn)  satisfying 

’  <  =  0  on  G0>  C  >  0  on  Rn  -  Gq 

« 

o  s  c  <  l. 

We  now  claim  that 


(3.13)  max(v.t+H(Dv.),v.-pC)  ^  0  on  Rn  x  {0} 

in  the  viscosity  sense,  which  means  that  If  ^  is  a  smooth  test  function  and 


(3.14) 


va-0  has  a  strict  local  minimum 
at  some  point  ( xQ , 0 )  €  Rn  x  (0>, 


then  either 
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(3. 15) 


v.(x-.O)  *  *i<(x  ) 


or  else 


(3. 16) 


bt  +  H(D*)  *  0  at  (xQ,0). 


Now  if  xn  €  Gn,  then  (3.15)  is  clearly  true.  Otherwise  suppose 


Xg  €  R  -  Gg  and 


v,(Xg,0)  <  MC(Xg)  <  09. 


We  repeat  now  the  argument  from  Steps  1-2,  noting  in  particular  that  since 


v^r(x,0)  =  +  oo 


for  all  x  near  Xg, 


the  points 


(xCr,tc') 


above  lie  in  R  x  (0,oo).  As  such  the  maximum  principle 


argument  leading  to  inequality  (3.12)  is  valid,  and  the  rest  of  the  proof 
proceeds  as  before,  yielding  at  last  the  inequality  (3.16). 


Since 


v®(x,0) 


as  e 


(k=l, • • • , m) 


if  Xg  €  Gg,  we  have 


To  see  that 


v„  =  0  on  G_.  x  (0>, 


(3.22) 


v.  ■  +  oo  on  (F  -Gn)  x  {0>, 


choose  any  point  xQ  e  R  -  GQ  and  suppose  instead 


(3. 17) 


v.(x,0)  <  oo  . 


Fix  8  >  0  and  then  define 


5  lx"x0l 

Y  ( x ,  t )  *  -  - - - At, 


for  A  =  A($)  to  be  selected  below.  Since  v#  is  lower  semlcontinuous, 
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(3.18)  v„  -  has  a  minimum  at  a  point  (xg,tg)  €  Rnx[0,eo). 


Then 


(3. 19) 


lvx0l2  IVxol2 

-  5  v*(*5.tg)  +  -  +  XtA  5  V,(xn,0)  <  00. 


Now  if  tg  >  0, 


whence 


+  H(D*5)  it  0  at  (xg, tg) ; 


(3.20) 


-A  + 


r  2(xa“xo)'i 

H[  -  —V2-)  2  0. 


a  contraction  for  A  =  A(S)  sufficiently  large.  Thus  tg  =  0.  If 

v*(xQ. 0)  <  mC(xq), 

then  (3.19)  Implies 

v.(xg,0)  <  MC(xg) 


for  small  enough  5,  and  so  according  to  (3.13)  we  once  more  would  obtain 
(3.20).  Thus 

v*(x0»0)  *  M<(x0). 

But  since  <(xQ)  >  0  and  p  >  0  is  arbitrary,  (3.17)  cannot  be  true.  □ 


Following  next  is  the  analogue  of  Proposition  3.1,  with  v  in  place  of 

v.- 

Proposition  3.2.  We  have 
(i) 

l 

(3.21)  min  (v*  +  H(Dv*),v*)  50  in  Rn  x  (0,«) 

in  the  viscosity  sense,  and 


» 
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(3.22) 


0  on  Gq  x  {0} 


+«  on  (Rn-GQ)  x  {0>, 


Proof  1.  Since  v  i  0,  to  establish  (3.21)  we  must  show 


•  *  • 

+  H(Dv  )  <  0  on  the  set  (v  >  0} , 


in  the  viscosity  sense.  So  select  any  smooth  test  function  <p  and  suppose 


(3.23) 


v  ~4>  has  a  strict  lc 
(xQ,  tQ)  €  IRn  X  (0,  ») , 


local  maximum  at  some  point 


(3.24) 


v  (xQ,t0)  >  0. 


We  need  to  show 


(3.25) 


<Pt  +  H(D*)  jS  0  at  (x0,tQ). 


Let  (^j,,,,,0jc)  be  a  positive  eigenvector  of  A(D^(xQ, tQ) ) ,  correspond¬ 
ing  to  the  principal  eigenvalue  A°( A(D^(xQ, tQ) )  =  H(D#(x0> tQ) ) .  Then 


(3.26) 


v  (x,t)  =  11m  sup  max 

e-)0  ISISm 
(y,s)-»(x,  t) 


s)  +  e  log 


Combining  (3.23)  and  (3.26),  we  deduce  that  there  exists  an  index 

G  G 

ke  {!,•••,*},  a  sequence  er  —*  0,  and  points  (x  r,t  r)  such  that 


(3.27)  v£rJxe,\tCrJ  +  er  lo6  ■  »x  [vir [xCr« tCp]  +  cr  lo8  “*  v 


(3.28)  v®p|xCrtterj  +  cr  *08  0k  “  0  has  a  local  maximum  at  £xC|\terj, 
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and 


(3.29) 


[xcMCr] 


(Xg.tg)  as  r  — >  oo. 


Utilizing  the  maximum  principle,  we  deduce  from  (2.6)^  the  inequality 


(3.30) 


f  uCr 

♦t  -  W*  *  dkl^|C  S  - 


at  the  point 


(*Cr •*'']• 


We  must  now  study  the  limiting  behavior  of  the  term 


on  the  right  hand  side  as  r  — >  oo. 


2.  We  assert 


(3.31) 


(1  =  1 ,  • • • , m) 


as  r  — )  oo.  To  see  this,  note  that 


exp  I- 


v  r  x  r  t 
vk  x  ,t 


owing  to  (3.24)  and  (3.27).  This  establishes  (3.31)  for  1  =  k.  Now  suppose 
additionally  that  for  some  index  1  *  k, 


lim  sup 


u^(xcMc-] 


a  >  0. 


Then  passing  if  necessary  through  an  appropriate  subsequence 


KL c  KL  •  * nnd 


fk[uc-(xc-,tG»)j  fk(  ,...,a,-**,0, ••.,), 


a  >  0  occupying  the  l^*1  argument  and  0  occupying  the  kfc^  argument  of  f. 
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Using  hypothesis  (F2),  we  obtain 


f.  xfc\t* 
11m  — — - 

G s  Go  , Gi 

S-XO  uk"  X  ,t 


=  +  CO, 


a  conclusion  at  variance  with  (3.30).  Thus  statement  (3.31)  is  true. 
And  Indeed  the  exact  reasoning  above  proves  additionally 


(3.32) 


“Cr 


=  0(1) 


as  r  — >oo,  the  functions  evaluated  at  the  point 


(*‘v.tCr]- 


We  further 


refine  this  observation  by  claiming 


(3. 33) 


fk 


C  r 

ckiV 


+  0(1) 


as  r  — »  oo.  To  see  this  we  observe  that 


(3.34) 


fk(u)  =  fk(0)  +  cklUl  +  °( lu|Z) 


=  cklUj  +  0(  |u|2)  (u  e  lRm). 


Now  (3.31)  yields 


m  cr 

->  u,r 


I  I  V'  1 

—  =o(1)  L  -i 

i «  ■  ii  * 


S  0(1) 


cr 

cklul 


1=1  V 


as  r  — »  oo.  Consequently  (3.34)  Implies 


fk  uGr  CklUir 

JL1 -  ;>  (i-o(l))  — 


as  r  — »  oo.  In  view  of  (3.32)  then 
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(3.35) 


=  0(1) 


£r 

CklUl 


as  r  — >  ».  Since 


Mu  1  ckiuir 

Cr  Er 

V  V 


according  to  hypothesis  (F5),  claim  (3.33)  is  proved. 


Finally  note  that  (3.35)  Implies 


(3.45) 


=  0(1) 


( 1=1 , • • • , m) 


as  r  — >oo. 


3.  Owing  to  (3.27)  and  (3.33),  inequality  (3.30)  yields 


*t  *dk 


as  c  — »  oo 

r 


at  the  point  (x0>tQ).  We  now  conclude  as  in  the  proof  of  Proposition  3 . 1 ( i ) 
4.  Next  we  verify  assertion  (il).  We  first  claim 


(3.37) 


«  *  * 

min  (vt  +  H(Dv  ),v  )  S  0  on  x  {0} 


in  the  viscosity  sense,  which  means  that  if  0  is  a  smooth  test  function  and 

* 

v  -<p  has  a  strict  local  maximum  at  some  point 

« 

(x.,0)  e  Gn  x  (0> 

0  0 

then  either 


(3.38) 


v  ( Xq , 0 )  =  0 


or  else 


(3.39) 


<f>t  +  H(D*)  S  0  at  (x0,0). 
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Now  if  (3.38)  is  false,  then  we  repeat  the  argument  from  Steps  1-3,  above 
noting  that  since 

v£(x,0)  =  0  for  all  x  near  xQ1 

the  points  |xCr,tCrj  lie  in  Rn  x  (0,oo).  Consequently  the  maximum  principle 
arguments  employed  above  lead  us  to  (3.39). 


5.  Now  observe  that 

v  =  +  oo  on  (I^-Gq)  x  {0>, 

£ 

since  v^  =  +  m  on  that  set.  Suppose  then  that  x^  6  Gq,  but 
(3.40)  v*(xQ, 0)  >  0. 


Fix  5  >  0  and  write 


x-x„ 


S 


+  At 


* 

for  A  =  A(5)  to  be  chosen.  Since  v  is  upper  semi continuous  and  is 
bounded  near  x^, 


(3.41) 


*  6 

v  has  a  local  maximum  at  a  point 

•  (xg. tfi)  €  Gq  x  [0, ») 


for  each  sufficiently  small  5  >  0,  with 


If  t_  >  0,  then 

O 


whence 


x5  XQ  as  3  -»  0. 


4>*  +  H(D*5)  SO  at  (x5,t5); 


(3.42) 


.  „(fal)  *  0, 
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a  contradiction  for  X  =  A(6)  large  enough.  Thus  t.  *  0. 

O 

Now 

0  <  v-[Vo)  S  v-(Xj,o]  -  . 

owing  to  (3.37)  we  once  again  reach  the  contradiction  (3.42).  Hence  (3.41) 
is  untenable,  so  that 

* 

v  =  0  on  Gq, 

as  required.  a 

Conclusion  of  the  proof  of  Theorem  1. 

In  light  of  Propositions  3.1  and  3.2  we  may  Invoke  the  uniqueness 
theorems  developed  in  the  Appendix,  §4,  to  find 

v  =  v,  =  I  in  Rn  x  (0,co), 

where  I  is  the  unique  viscosity  solution  of  the  Hamilton- Jacobi  variational 
inequality 


min(It+H(DI), I)  =  0 

in 

IRn  x  ( 0 ,  co) 

(3.43) 

A 

I  =  0 

on 

G0  x  <0> 

I  =  +  00 

on 

(Rn-G0)x{0> 

Additional ly, 

(3.44)  v£  — »  I  uniformly  on  compact  subsets  of  Rnx(0,«),  k=l,***,m. 

Now  according  to  f5  in  [6],  we  have 

(3.45)  I  =  max  (J.0), 

where  J  is  the  unique  viscosity  solution  of  the  Hamilton- Jacobi  equation 
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Jt 

+  H(DJ)  =  0 

in 

SRn  x  (0,oo) 

(3.46) 

J 

=  0 

on 

GQ  x  {0> 

J 

=  +  00 

on 

(Rn-GQ)x{0> 

In  addition  we  have  the  representation  formula 

(3.47)  J ( x , t )  a  inf  j  J  L(z(s) )ds|z(0)  €  Gg,  z(t)  =  xj, 

the  inflmum  taken  over  all  absolutely  continuous  functions 

z:  [0,  t }  -»  !Rn. 

Now  since  (2.5)  and  (3.44)  imply 

-  -I+o(l) 

u£  =  e  e  =  e  C  (k=l, • • • , m) , 

we  see  that 

(3.48)  u£  — »  0  uniformly  on  compact  subsets  of  (I>0>  =  (J>0> 

els  e  — »  0,  for  k=l, ••*,!»>. 

We  must  now  show 

(3.49)  lim  inf  min  u  >0  uniformly  on  compact  subset  of  {J<0>. 

e-)0  ISkSm 

So  fix  any  point  (Xg.tg)  s  {J<0>.  Then  1=0  near  (Xg.tg). 

♦(x.t)  *  h-x0|2  +  |t-tQ|2 
Owing  to  (3.44)  we  see  that  for  each  k=l,***,m. 

G  I  G  G 1 

(3.50)  vk  -  0  has  a  maximum  at  a  point  x^. tjJ  » 
with 

(3.51)  (xk’t^)  — >  813  c  —* *  °- 


Define 
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Applying  the  maximum  principle,  we  find  using  (2.6)e  that 


(3.52) 


and  so 


(3.53) 


as  e  — »  0. 


2  fk  U 

o(l)  =  *t-edkA*  +  dk|D*r  S  - 

Uk 


fk(uc)  S  o(l)  u£  at  the  point  • 


Now  there  exists  a  constant  a  >  0  such  that 


(3.54) 


fk(u)  >  c^u^  -  a|u|^  for  all  u  e  Rm 


Let  us  suppose  first  that 


(3.55) 


ui[xk-tk) 


for  each  1=1 , • • • , m, l*k. 


Then  from  (3.53)  and  (3.54)  we  deduce 


o(1K 2  fk,uC)  *  Z  cki“I *  °kk\ 

l*k 

(Xk’^k]:  w^ence  (3.55)  implies 
„ ,  e f  e  .  el  .  ckk 

d)  \  [xk’  tkJ  *  2a 


(3.56) 


for  e  small  enough.  Should  (3.55)  fail,  then 


ef  e  .el  ckl  ,, 

ullxk’tkJ  >  2a~  for 


some  1  *  k. 


But  then  owing  to  hypothesis  (F2) 


at  the  point 


(xk*  lk)  * 


f  (-Vi,0iVi-",)  *  0 


for  some  positive  constant  p.  Thus  at 


(<•‘3 
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have 


fk(ue)  u  +  fk(ue)  -  fk(***uJ_1.0. V!--’*) 

4 8 *  ckk\  -  '[M2  *  z  v?] 

1*1 


for  some  y  >  0.  Consequently  (3.53)  implies 

[uk]  *  8  *  °(Kl)  at  (xk’tk)’ 

and  so 

uk  [<■  lk) 2  6  >  0 

for  sufficiently  small  5  >  0.  But  since 


we  have 


and  so 


Thus 


4.  Appendix: 


vk(^’tk]  2  Vk(xO'to): 

'iH-1*}  S  Uk(X0>to)- 

lim  inf  u£  ^x0> Iq]  *  s  >  °- 
Identification  of  the  action  function 


a 


We  outline  in  this  section  a  proof  that  the  functions  vm  and  v 
Introduced  in  f3  agree  and  equal  I,  the  unique  viscosity  solution  of 


(4. 1) 


min(It+H(DI),I)  =  0 

in 

Rn  x  (0,«) 

I  =  0 

on 

GQ  x  (0> 

8 

♦ 

II 

on 

(Rn-G0)x{0> 
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First  we  show 


(4.2)  v0  2  I  in  Fn  x  (0,eo). 

For  this  choose  any  constant  y  >  0  and  any  function  <  €  C®(Rn)  satisfying 
'  <  a  0  on  GQ,  C  >  0  on  Rn  -  GQ 

* 

_  0  iS  <  £  1. 


Consider  now  the  auxiliary  problem 


(4.3) 


f  min(I  . +H(DI  ) , I  )  =  0 
1  yt  y  y 


.  V  = 


in  Rn  x  (0, ») 
on  Rn  x  (0,oo)  , 


which  has  a  unique,  continuous  solution  according  to  [6].  Furthermore, 
according  to  assertion  (ii)  in  Proposition  3.1  we  have 

(4.4)  vm  2  I  on  Rn  x  (0). 


Finally,  vm  is  lower  semicont inuous,  and 

min(v.t  +  H(Dv.),v„)  20  in  Rn  x  (0,») 


in  the  viscosity  sense.  Consequently  a  comparison  argument 
following  [3]  and  [6]  implies 

v,  2  I  in  Rn  x  (0,  co) . 

Letting  y  tend  to  infinity  we  have 

I  — »  I  in  Rn  x  (0,  co) ; 

and  so  (4.2)  follows. 

Next  we  assert  that 

(4.5)  v  £  I  in  Rn  x  ( 0, co) . 
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To  prove  this  define  for  each  small  5  >  0  the  smooth  set 
Ga  ■  {xQ  €  Gq |dlst (x, Rn-GQ)  >  5>. 

Fix  8  >  0  and  write 

(4.6)  A  *  sup  v  (x,«r)  (<r  >  0). 

'  X€G5 

Since  v  is  upper  semi continuous  and 

v  =0  on  Gq  x  {0}  , 

according  to  assertion  (ii)  of  Proposition  3.2,  we  see  that 


(4.7) 


lim  A  =  0 
<r-i 0  a 


for  each  fixed  8  >  0.  Choose  some  small  <r  >  0  and  consider  the  problem 


(4.8) 


minfl^,<r+H(PI5,<r),  I5,<r)  =  0  in  Rn  x  (<r,eo) 


5,<r  _ 


+00 


In  view  of  Proposition  3.2  we  have 
v  <  I 


on  Gfi  x  {<r> 
on  (Rn-Ga)x{o-> 


on  Rn  x  {o-} 


Since  additionally  v  is  upper  semlcontinuous  and 


*  *  * 

mln(v  +  H(Dv  ),v  )  SO 


in  the  viscosity  sense,  we  have 

*  *  T5-0* 

v  <  I 


in  Rn  x  (a,  eo ) 


in  Rn  x  (o-, oo) , 


Let  o'  — >  0  and  recall  (4.7)  to  discover 

(4.9)  v*  S  I5  on  Rn  x  (0,»), 


when 


is  the  unique  viscosity  solution  of 
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f  min fl^  +H(DI6),!6)  =  0 


Now  at  last  send  S  — »  0:  since 


in  SRn  x  (0,oo) 
on  Gg  x  <  0> 
on  (Rn-G5)x{0>  . 


in  Rn  x  ( 0 ,  oo ) , 


we  arrive  at  (4.5). 

Combining  (4.2)  and  (4.5)  we  have 

v  £  I  £  v„ 


in  Rn  x  (0, «) . 


But  since  the  definitions  imply  obviously  that 

v.  <  v*  in  Rn  x  (O.oo), 

we  have 

v*  =  v.  =  I  in  Rn  x  (0, oo) .  □ 


Another  approach  to  obtain  the  above  is  to  modify  (in  a  more  or  less 
straight  forward  way)  the  results  of  H. G.  Crandall,  P.-L. Lions  and  P.E. 
Souganidis  [4]  concerning  maximal  solutions. 
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